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Abstract 

Using the renormalized quasiparticle random phase approximation (RQRPA), 
we calculate the light neutrino mass mediated mode of neutrinoless double 
beta decay (Oz^/?/3-decay) of ^^Ge, ^^Mo, ^"^^Te and ^^^Te. Our results indi- 
cate that the simple quasiboson approximation is not good enough to study 
the Oz^/3/3-decay, because its solutions collapse for physical values of gpp. We 
find that extension of the Hilbert space and inclusion of the Pauli Principle in 
the QRPA with proton-neutron pairing, allows us to extend our calculations 
beyond the point of collapse, for physical values of the nuclear force strength. 
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As a consequence one might be able to extract more accurate values on the 
effective neutrino mass by using the best available experimental limits on the 
half-life of Oi^/3/3-decay. 
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Neutrinoless double beta decay (Oz//9/3), which involves the emission of two electrons and 
no neutrinos, requires the neutrino to be a Majorana particle with non-zero mass. This 
process violates the lepton number conservation and occurs in some theories beyond the 
standard model (see e.g. for reviews). The Oi//3/5-decay is not yet observed. The 

experimental lower limits on the half-lives of some nuclei which undergo double beta decay, 
provide us with the most stringent limits on the effective mass of the electron neutrinos 
and the parameters of the right-handed currents after evaluating the corresponding nuclear 
matrix elements. 

The calculation of the Oi//5/5-decay matrix elements is based on some approximation 
schemes, from which the Quasiparticle Random Phase Approximation (QRPA) is the most 
prominent one 0-0, despite an uncertainty in the calculation of the nuclear matrix element, 
which is related to its sensitivity on the renormalization of the particle-particle component 
of the residual interaction. Thus in the framework of proton-neutron QRPA (pn-QRPA) 
using the zero range delta-force Engel et al. have found a strong suppression of the Qv(3(3- 
decay matrix elements by introducing short-range correlations, whereas Tomoda et al. 
and Muto et al. |]^ by using a more realistic effective NN interaction arrived at the opposite 
conclusion. 

Recently new approaches for the calculation of the nuclear many body Green function 
have been proposed, which are supposed to offer more reliable results than that of pn-QRPA. 
The pn-QRPA has been extended to include proton-neutron pairing (full-QRPA) 0. It has 
been found that proton-neutron pairing influences the Oz//5/5-decay rates significantly HlO 



Toivanen and Suhonen [|TT[], Schwieger, Simkovic and Faessler [|1^ and F. Krmpotic et al. 



131 have studied the effect of the violation of the Pauli principle in the correlated ground 



state of the pn-QRPA introducing the renormalized pn-QRPA (pn-RQRPA). Schwieger et 



al. have done calculations of the two-neutrino double beta decay matrix elements using 
this renormalized pn-QRPA but also including proton- neutron pairing (full-RQRPA). They 
have shown that the inclusion of both p-n pairing and an improved treatment of the 



Pauli principle removes the difficulties arising from the strong dependence of the matrix 



element on the particle-particle strength Qpp in the standard pn-QRPA. 

The Ou (3 jS-decaj allows to determine an upper limit of the effective Majorana electron 
neutrino mass < rrii, > [|1|- [Q. For this process the most stringent experimental lower limit 
of the life-time for the Oi//3/5-decay of "^^Ge favours especially this nucleus for extracting an 
upper limit of < m,^ >. The purpose of this article is to study in detail the effects of the 
Pauli principle, the proton-neutron pairing and the truncation of the nuclear Hamiltonian in 
the evaluation of the nuclear matrix element governing the Oz//3/3-decay of '^^Ge. In addition 
the fuU-RQRPA will also be used to calculate 0u(3f3-decay of ^°°Mo, ^^^Te and ^^^Te. The 
stability of the nuclear matrix elements in respect to the renormalization of the particle- 
particle force will be discussed. 

The fuU-RQRPA, which describes the excited states of the even-even nucleus, has been 
studied in Ref. [1^. Therefore, here we shall present only the formulae relevant to this work. 



In the notation of Ref. |]T2[, the half life of the Oi//?/3-decay in the case of the neutrino 



mass mechanism is given by, 

[Ty,]-^ = Goi(MlJ^(^^)' (1) 

with the effective neutrino mass 

< >= ^ |f/,/m,e'"% (2) 
j 

where exp(iaj) is the CP eigenvalue of the neutrino mass eigenstate {uj >, Uei is the element 
of the unitary neutrino mixing matrix and Trig is the mass of the electron. Gqi is the 
integrated kinematical factor for the Oj^ —>■ Oj transition [^p!0|. The nuclear matrix 
element M2!f„^^ consists of Fermi and Gamow- Teller contributions 



mas 



<l. = M°^^-[^yMr (3) 



where 



■ 'H{n,) ), (4) 

M^T I ^ CTi • (T2 



< Oi2 > = E {-Y+^>''+'+^{2J + 1) 



jk jl J 

I ^ 

X < pk,pk';J\f{n2)r+40i2f{n2)\nl,nl';J> 
X < 0/ II [Cpk'^^i'^J II ^^"^/ >< J^m/I J'^mi >< J'^m, || [cj.c„i]j || 0+ > . (5) 

The short-range correlations between the two interacting nucleons are taken into account 
by a correlation function f{r 12) |§|^- The neutrino-potential H{ri2) in the case of light 
neutrino exchange takes the form 

rrf \ _ ^ f°° sin(gr) 1 

Here, R = tqA^^^ is the nuclear radius (rg = 1.1 fm) and the parameter A of the dipole 



shape nucleon form factor is chosen to be 0.85 GeV ||8|JT0|]. In the fuU-RQRPA for the 



matrix elements of the one-body transition densities we obtain 



< rm, II [c;,Cru]j II 0+ > = V2JT1 E rnif^k^ul) [utl.vlllX ^^ikj, r) 

/X,!/ = l,2 



<Q^\\[4,rcnv]j\\rmf> =V2jTi E ni{^ik\ui')[v\!,l^u\!lx';i{k\i\r) 

^J,,l^=l,2 



with m(/ia, ub) = ^+1"^^^^/^ . We note that the X^^{k, I, J"") and F^^(A;, /, J"") amplitudes 
are calculated by the renormalized QRPA equation only for the configurations fia < ub { 
i.e., n = u and the orbitals are ordered a < b and /i = 1, z/ = 2 and the orbitals are not 
ordered) ||12|. For different configurations 'X^^ik^l, J'^) and Y^j^{k,l, J'^) in Eqs. (|^) and 
(1) are given by following the prescription in Eqs. (65) and (66) of Ref. P|. The index i 
(f) indicates that the quasiparticles and the excited states of the nucleus are defined with 
respect to the initial (final) nuclear ground state |0j^ > (|0j >). The overlap between two 
intermediate nuclear states belonging to two different sets is given in Ref. ||12[. 
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We note that for Vk^iuj^ = 1 the expressions (J^) and (H) are just the one body transition 
densities in the fuU-QRPA |T^. If in addition U2p = V2p = Win = ^in = (i.e. there is no 
proton-neutron pairing), Eqs. (0) and (|[) reduce to the expressions of the pn-QRPA • 

We apply the pn-QRPA, fuU-QRPA (with proton- neutron pairing), pn-RQRPA and fuU- 
RQRPA methods to the 0z//3/?-decay of "^^Ge to study the effects of the proton-neutron 
pairing and the Pauh-principle violation on the nuclear matrix element M^^^^g. In previous 
calculations that model space has been used, which described best the beta strength distri- 
butions for the initial and final nuclei. However, there is a substantial difference between the 
beta strength matrix element and the double beta decay matrix element. Unlike the case 
of the double beta decay, which is a second order process, meson exchange currents do not 
play an important role in the calculation of the single beta strength distribution. The meson 
exchange currents are incorporated in our calculation through the G-matrix elements of the 
nuclear Hamiltonian. We suppose that they can be taken into account in a proper way only 
if the inclusion of Pauli principle is considered, and a large enough model space is assumed. 
In addition we wish to stay as closely as possible to the Brueckner reaction matrix of the 
Bonn potential. A smaller model space requires a more drastic renormalization of the force. 
Since one does not know a reliable microscopic theory for the renormalization as a function 
of the model space, a smaller model space means more free parameters. Thus we enlarge 
the model space of the nuclear Hamiltonian. The model spaces used are the following: 
9 — level model space : The full 3huj and Ahoo major oscillator shells. This model space has 
been used by Tomoda et al. ||^ and by Muto et al. 0. 

12 — level model space : It consists of the full 2 — Ahu major oscillator shells. 

21 — level model space : This model space extends over the — 5huj major oscillator shells. 

The single particle energies have been calculated with a Coulomb-corrected Woods-Saxon 
potential. As a realistic two-body interaction we use the nuclear matter G-matrix calculated 
from the Bonn one-boson-exchange potential |T^. We introduce the parameterization daG 



(a=pp, nn and pn), gphG and gppG for the pairing, particle-hole and particle-particle interac- 
tion, respectively. The quasiparticle energies and amplitudes have been found by solving the 
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HFB-equation with proton- neutron pairing . For a model space larger than two oscillator 
major shells we neglect the mixing of different "n" but the same "Ijm" orbitals. We suppose 
that the Woods-Saxon potential is already a good potential and therefore shell mixing is not 
appreciably affecting the double beta decay. The renormalization of the pairing interaction 
has been determined to fit the empirical proton, neutron and proton- neutron pairing gaps 
according to Ref. The renormalization parameter da together with the experimental 
proton ( Ap^'^ ), neutron ( /S.^^'^ ) and proton - neutron ( ) pairing gaps for ^^Ge and 
'^^Se are listed in Table |. We see that dpp and dnn values are close to unity. A dpn value 
higher than unity is the price paid for the spherical symmetry of the model which excludes 
the treatment of the T=0 pairing. The J=0 T=0 pairs can be treated in a BCS or even 
HFB approach only due to deformation. The T=0 pairing is effectively taken into account 
by the renormalization of the T=l J=0 n-p interaction leading to a higher value of dpn. We 
note that in the framework of the HFB method we have found that the importance of the 
proton-neutron pairing for the ground state properties of nuclei is decreasing significantly 
as isospin increases P]. This fact is in agreement with the recent sophisticated study of J. 



Engel et al. [1T3 



Further, we fixed Qph = 0.8 as in our previous calculations and only the particle- 
particle interaction strength Qpp is considered as a variable. Fig. 1 (a) shows the calculated 
nuclear matrix element M^^^^ in the framework of the pn-QRPA for the Ou/SjS of ^^Ge 
as function of gpp. It is worthwhile to notice that by increasing the model space M^^^^ 
becomes extremely sensitive to the Qpp and even within the physically acceptable region of 
Qpp collapses and crosses zero. This behaviour has its origin in the contribution of the J'^ = 1"*" 
intermediate states of the odd-odd nucleus to the many body Green function M^'^^g, which 
becomes too large because of the generation of too many ground state correlations. This 
feature has not been found in the previous calculations as the model spaces used there, were 
too small. We surmise that the different suppressions of the Oz//3/3-decay matrix elements 
found by Tomoda et al. [^] and Muto et al. [0 on one side and Engel et al. on the other 
side could have their origin not only in the different interactions but also in the different 



model spaces used. Unfortunately, the model space of Ref. is not specified to draw a 
definite conclusion. F. Krmpotic and S. Sharma [0] have calculated OujSP-decaj of ''^Ge 



by using eleven dimensional model space (The full 3huj and 4huj major oscillator shells plus 
0/2,9/2 and 0/19/2 orbitals.), zero range 5-force interaction and different treatment of the two- 
nucleon correlation function. Their results show also a strong suppression of the M^^^^ with 
increasing strength of the particle-particle force. To our knowledge, we are the first to study 
in details the extension of the Hilbert space for the Oi>(3(3-decaj of ^^Ge and showing that its 
solution collapse for physical value of gpp in the QRPA. Therefore, pn-QRPA method does 
not allow us to make definite Oz//3/3-decay rate predictions. 

In Fig. 1 (b) we show the results obtained in the framework of the full-QRPA, which take 
into account proton-neutron pairing [|10[. The M^^^^ has been calculated till the collapse 



of the full-QRPA, which occurs for the J'^ = 0^ or channel. The early collapse of the 
full-QRPA for the 21-level model space does not allow us to draw conclusions inside the 
expected physical range 0.8 < Qpp < 1.2 within this approximation. 

The nuclear matrix element M^^^^ obtained within the pn-RQRPA is shown in Fig. 1 
(c). From the comparison with the Fig. 1 (a) it follows that the inclusion of ground state 
correlations beyond the QRPA in the calculation of M^^^^ removes the difficulties associated 
with the strong dependence of M^'^g^ on the particle-particle strength. It means that the 
Pauli principle plays a prominent role in the evaluation of the nuclear matrix elements 
governing the second order processes. We conclude, that the rather big difference between 
the results of 9-level and 12-level model spaces and the small difference between the results 
of 12-level and 21-level model spaces favours the 12-level model space for the study of double 
beta decay of ''^Ge. 

In Fig. 1 (d) we present the results with the full-RQRPA method. This method takes into 
account both proton-neutron pairing and the Pauli exclusion principle and does not show a 
breakdown for any realistic interaction strength. A comparison with Fig. 1 (c) shows that 
proton-neutron pairing decreases the values of M^^^^ for the small model space (9-levels) 
and increases them for bigger, more realistic model spaces. Thus the full-RQRPA allows an 
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extension of the calculations beyond the point of collapse and supports more stable solutions 
for calculating the M^^^^ , with respect to the renormalization of the particle-particle force. 

In order to gain more confidence in the fuU-RQRPA method, we have also performed 
calculations for Oi//3/5-decay of ^°°Mo, ^^^Te and ^^°Te. As noticed already in Ref. |]T0|, the 



evaluation of M^^^^ is very sensitive to the truncation of the model space in the case of 
^^^Mo. We have found that the model spaces used in the previous calculations were not 
sufficient to stabilize the results and therefore we have performed calculations with the 21- 
level model space defined above. For the Te isotopes we have introduced a model space, 
which comprises 20 levels: The full 2-5 hu major oscillator shells and 0211/2,13/2 subshells. 
The calculated matrix elements M^^^^ are shown in Fig. 2. For all nuclei studied M^^^^ 
is only weakly dependent on the strength of particle-particle interaction. This allows us to 
have more confidence in the limits for the effective electron neutrino mass < >extracted 
by using the best presently available experimental limits on the half-lives of the Oi//3/?-decays. 
We obtain the following limits: 
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We see, that by increasing significantly the model space for ^°°Mo the value of M^^^^ is much 
bigger as in our previous calculations [|10|. The rather big energy release for the Oz//3/3-decay 



of ^°°Mo and the large value of the matrix element M^'^^g favour this nucleus for further 
experimental study. 

The above results suggest, that in the calculation of the nuclear matrix element M^^^^, 
it is necessary to take into account the Pauli principle and to consider a large enough model 
space in addition to pn-pairing. As it was noticed already by J. Hirsch, P. Hess and O. 
Civitarese |PT| and by F. Krmpotic, A. Mariano, E.J.V. Passos, A.F.R. de Toledo Piza and 



T.T.S. Kuo H, the price which is paid for taking into account the ground state correlation 
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beyond the QRPA, is the violation of the Ikeda Sum Rule. In the framework of the full- 
RQRPA we have 

S--S+ = ^ I < A; II a II / > |2 

kfilu 

^ (^/c/^n '^kfip ~^ '^kfipUk^n'^lupUlun UkfMpVk^nUlupVlun) ■> (9) 

which yields 3(N-Z) only when D]^^iyi+ = 1 (quasiboson approximation). If the usual 
HFB constraint on the particle number is used and if the model space contains both states 
j+i/2,-1/2 = /±l/2 of the spin-orbit splitting, the QRPA fulfills the Ikeda Sum Rule indepen- 
dent of the chosen model space, which is not a realistic feature of the model. In the RQRPA 
we necessarily face its violation because of D 1.^1^1+ 7^ 1. By using the full- RQRPA for "^^Ge 
and ^^Se we have found — to be 0.8 — 0.9 x 3(A^ — Z) in the physical acceptable 
region of the parameter Qpp. It is supposed that the omission of the scattering terms in the 



operators for and /3~ transitions, could be the reason for this small mismatch ||T3|,^. It 
is worthwhile mentioning that the Ikeda Sum Rule is not satisfied even if the BSC equations 
are solved within the RQRPA with the condition that the average particle number in the 



correlated ground state is conserved |T^. Nevertheless we expect that a small violation of 
the Ikeda Sum Rule wont alter our main results. 

In summary, we have studied the effects of Pauli principle violation and proton-neutron 
pairing on the calculation of the Oz//5/3-decay matrix element. We have shown that the 
renormalised QRPA allows extending the calculations beyond the collapse and thus that the 
simple quasiboson approximation is not sufficiently accurate to calculate reliable the nuclear 
many-body Green function governing the Oz//?/5-decay process. This becomes evident, if a 
larger, more realistic, Hilbert space is considered. We have found that the Oz//?/5-decay 
matrix element calculated via the fuU-RQRPA, which includes the Pauli effect of fermion 
pairs, avoids the collapse for physical values of the nuclear force, and therefore we can argue 
that we have some evidence that the Pauli effect violation inherent in the traditional QRPA 
affects the Oz//3/3-decay matrix element significantly. As a consequence, to our opinion, the 
results presented in the text for the effective neutrino mass for Qpp = 1 are more accurate. 
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Note added. After completing this work, we recieved a preprints by J. Hirsch et al. 



22 1 and by J. Engel et al. p3| who studied the validity of the renormalized QRPA within 



a schematic exactly solvable models, which are not intended to reproduce actual nuclear 
properties. They confirmed that the renormalized QRPA offers advantages over the QRPA. 
However, they founded some discrepancies between the exact and the RQRPA solutions 
after the point of collapse of the QRPA. It is questioned whether these discrepancies should 
hold in more realistic calculations. If they would, it just supports our conclusion that we 
need a theory in which the Pauli principle is explicitly incorporated. In the renormalized 
QRPA the Pauli principle is considered in an approximate way. We note that there is no 
exactly solvable realistic model. 
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TABLES 

TABLE I. Experimental proton ( A^^^ ) , neutron ( A^^^ ) and proton - neutron ( ) pairing 
gaps and renormalization constants of the proton - proton ( dpp ) , neutron - neutron ( dnn ) and 



proton - 


neutron ( dpn ) pairing 


interactions for studied nuclei. 






Nucleus 


( A'^^P A<^^P A'^^P 
( [MeV] ) 


model 
space 


7 

Upp 


dnn 


dpn 




(1.561, 1.535, 0.388) 


9 level 


1.063 


1.238 


2.093 






12 level 


0.988 


1.150 


1.777 






r\ -1 1 1 

21 level 


0.899 


1.028 


1.506 




(1.751, 1.710, 0.459) 


9 level 


1.135 


1.255 


1.678 






12 level 


1.027 


1.182 


1.524 






21 level 


0.934 


1.059 


1.325 




(1.612, 1.358, 0.635) 


21 level 


0.980 


0.923 


1.766 




(1.548, 1.296, 0.277) 


21 level 


1.002 


0.945 


1.568 




(1.127, 1.177, 0.149) 


20 level 


0.873 


0.942 


1.780 




(1.043, 1.180, 0.090) 


20 level 


0.835 


0.945 


1.87 


IfXeu 


(1.307, 1.266, 0.199) 


20 level 


0.920 


0.972 


1.530 




(1.299, 1.243, 0.190) 


20 level 


0.914 


0.974 


1.614 
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FIGURES 

FIG. 1. The calculated nuclear matrix element M^^^^g for the Oz^/?/?-decay of ^^Ge as a function 
of the particle-particle interaction strength gpp. In (a), (b), (c) and (d) M^^^^ have been calcu- 
lated within the pn-QRPA, full-QRPA, pn-RQRPA and full-RQRPA, respectively. The dotted line 
corresponds to the 9-level model space, the dashed line to the 12-level model space and the solid 
line to 21-level model space, respectively. 

FIG. 2. The calculated nuclear matrix element M^^g^ for the Oz^/3/3-decay of ^^Ge within the 
full-RQRPA as function of particle-particle interaction constant gpp. In the case of ^^Ge and ^°°Mo 
(^^^Te and ^^"Te) the 21-level (20-level) model space has been used. 
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